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Abstract. 

We estimate sensitivity coefficients to variation of the electron-to-proton mass 
ratio n for microwave transitions in partly deuterated ammonia NII2D and ND2H. 
Because of the mixing between rotational and inversion degrees of freedom the 
coefficients Q^t strongly depend on the quantum numbers of the transition. This can 
be used for astrophysical search for possible variation of the constant fi. 
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1. Introduction 

At present discrete microwave spectra of molecules are used for astrophysical studies 
of possible variation of the fine structure constant a — e^/(^c), the electron-to-proton 
mass ratio /i — me/mp, and the nuclear g^- factor g^- Rapid progress in experiments with 
cold molecules can lead to new high precision laboratory tests of the possible variation 
of II (see, for example, review [1] and references therein). It was pointed out in [2] 
that inversion transitions in fully deuterated ammonia ^^NDs have high sensitivity to 
//-variation, = 5.6 (sensitivity coefficients Q are defined in section 2). The inversion 
transition in non-deuterated ammonia has a slightly smaller sensitivity to /i-variation, 
Qfj_ = 4.5 [3]. Molecular rotational lines have much smaller sensitivity, = 1.0. In 
astrophysics the observed frequency shifts are interpreted as Doppler shifts. Because 
of that, possible //-variation would lead to apparent velocity offset between ammonia 
inversion line and rotational molecular lines, originated from the same gas clouds. This 
fact was used in [3, 4, 5] to establish very stringent limits on /(-variation over cosmological 
timescale ~ 10^° years. 

Recently ammonia method was applied to dense prestellar molecular clouds in the 
Milky Way [6, 7, 8]. These observations provide a safe bound of a maximum velocity 
offset between ammonia and other molecules at the level of \AV\ < 28 m/s. This bound 
corresponds to | A/x/ /t| < 3 x 10~^, which is two orders of magnitude more sensitive than 
extragalactic constraints [3, 4, 5]. Taken at face value the measured AV shows positive 
shifts between the line centers of NH3 and other molecules and suggests a real offset, 
which would imply a Aji/ ji — (2.2 ± 0.4stat ± 0.3sys) x 10~^ [8]. These results can be 
relevant to the theories, which predict dependence of the fundamental constants on the 
local matter density [9]. 

One of the main possible sources of the systematic errors in such observations 
is the Doppler noise, i.e. stochastic velocity offsets between different species caused 
by difi^crcnt spacial distributions of molecules in the gas clouds [10, 11]. Because of 
that it is preferable to use lines with different sensitivity to variation of fundamental 
constants of the same species. Recently it was shown that sensitivity coefficients for 
A-doublet spectra of OH and CH molecules strongly depend on quantum numbers 
[12, 13]. In this context partly deuterated ammonia molecules NH2D and ND2H may 
be also interesting. Due to the broken symmetry, the rotational and inversion degrees 
of freedom for these molecules are strongly mixed. As we will show below, this leads to 
a significant variation of the sensitivity coefficients of different microwave transitions. 
Note that microwave spectra of NH2D and ND2H from the interstellar medium were 
recently detected [14, 15, 16, 17, 18, 19, 20, 21, 22]. 
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2. Sensitivity coefficients 

Let us define dimensionless sensitivity coefficients to the variation of fundamental 
constants so that: 

^ = C.^ + 0„^ + <2,^. (1) 

UJ a II Qn 

These coefficients Qi are most relevant in astrophysics, where lines are Doppler 
broadened and hnewidth T — ux SV/ c, where 6V is the velocity distribution width 

and c is the speed of light. The redshift of a given line is defined as Zi = uj\^h,il'^i — 1- 
Frequency shift (1) leads to the change in the apparent redshifts of individual lines. The 
difference in the redshifts of two lines is given by: 

^i^ = -AQ,:?^^-Ae„^-A(J,^»!l. (2) 

1 + z a IX 

where z is the average redshift of both lines and AQq, = Qa,i — Qa,j, etc. We can rewrite 
(2) in terms of the variation of a single parameter J-': 

= J-=a^'3>^%^««. (3) 



1 + z T 

The typical values of SV for extragalactic spectra is about few km/s. This determines 
the accuracy of the redshift measurements on the order of bz — 10~^ - 10~^, practically 
independent on the transition frequency. For gas clouds in the Milky Way the accuracy 

can be two orders of magnitude higher, bz = 10~^ ~ 10~^. In both cases the sensitivity 
of astrophysical spectra to variations of fundamental constants directly depends on AQj. 

In the optical range the sensitivity coefficients arc typically on the order of 10^^ 
- 10~^, while in the microwave and far infrared frequency regions Qi ~ 1. However, 
Equation (3) shows, that we need lines with different sensitivities. It is well known that 
for rotational transitions = 1.0, whereas for vibrational transitions — 0.5. For 
both of them \Qa\ <S 1 and \Qg\ <^ 1. Inversion transition in NH3 has Qfj, = 4.5. In the 
microwave region one can also observe hyperfine transitions {Qa — 2, Qf^ — 1, Qg — 1) 
and A-doublet transitions, where Qa and Q^ strongly depend on quantum numbers and 
can be very large [13]. This makes observations in microwave and far infrared wavelength 
regions potentially more sensitive to variations of fundamental constants, as compared 
to optical observations. Because of the lower sensitivity, systematic effects in the optical 
region are significantly larger [23]. 



2.1. Inversion transitions 

Sensitivity coefficients for inversion transitions in NH3 and ND3 were obtained in 
[2, 3]. For non-symmetric molecules NH2D and ND2H selection rules arc such that 
purely inversion transitions are not observable. Still, we will first estimate sensitivity 
coefficients for the inversion transition and then will proceed to the mixed inversion- 
rotation transitions. 
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Table 1. Sensitivity coefficients for the inversion transitions in different 
isotopologues of ammonia. 



Molecule 


Action 
S 


Q 

this work 


[2] 


[3] 


"NH3 


5.9 


4.4 




4.5 




6.0 


4.4 








6.5 


4.7 






"ND2H 


7.3 


5.1 






"ND3 


8.4 


5.7 




5.7 




8.5 


5.7 


5.6 





In the WKB approximation the inversion transition frequency is given by the 
expression [24]: 

oJinv = — , (4) 

TT 

where ujy is the vibrational frequency for the inversion mode and S is the action (we use 
atomic units h — rUe — e — 1). Vibrational frequency is inversely proportional to the 
square root of the reduced mass Mi and, hence, is proportional to /x^/^. The action. 



S = / ^J2Ml[U{x) - E]dx, (5) 

J —a 

to a first approximation, is proportional to (integration here goes between classical 
turning points, i.e. U{±a) — E). Therefore, we can rewrite (4) in the form [2]: 

f^inv = aijj}''^ exp (-02//"^^^) , (6) 

and find the respective sensitivity coefficient: 

Q^ = i(l + a2/i-i/2) =1(1 + 5) . (7) 

Deriving (7) we neglected that the energy E in (5) also depends on /j,, E — 
Umin + |<^v Taking this into account we get corrected expression [3]: 

Following [25] we can estimate AU = f/max ~ f^mm for ammonia to be approximately 
2020 cm~^. Now we can use experimental frequencies and cjjnv for different isotopic 
variants of ammonia to find S from (4) and estimate using (8). Results are presented 
in table 1. 

2.2. Mixed transitions 

For partly deuterated ammonia inversion levels have different ortho-para symmetry. 
Because of that inversion transitions can be observed only in combination with rotational 
transitions cUr. For such a mixed transition, 

U = LJ,:± Winv , (9) 
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and sensitivity coefficient is equal to: 

Qu — Qr,u, i Qinv.u -i \^^) 

UJ UJ 

where Qr,/* = 1 and Qinv,/* is given in table 1. 

Equation (10) shows that the sensitivity coefficient of the mixed transition is simply 
a weighted average of those of constituents. For ammonia all observed transitions 
have complex hyperfine structure. That means that in addition to two term in (9) 
there is third hyperfine term Whf. However, the hyperfine contribution is typically very 
small, |(^iif/^| ^ 1- Therefore, we can neglect hyperfine corrections to the sensitivity 
coefficients. Exceptions to this rule take place when there are accidental degeneracies 
between different levels and respective transition frequencies become comparable to 
hyperfine structure. This happens, for example, for some A-doublet transitions in NO 
and LiO molecules [13]. 



3. Molecules NH2D and ND2H 

Partly deuterated ammonia NH2D and ND2H are asymmetric tops. Their rotational 
levels are classified by the rotational quantum number J and two projection Ka and 
i^c; which correspond to two limiting cases of prolate and oblate symmetric tops. 
Thus, rotational levels are labeled as JKa,K^- Symmetric and asymmetric levels, which 
correspond to the inversion tunneling mode are labeled with indexes s and a respectively. 

Total inversion symmetry of the rotation- inversion levels is (—1)^'' for s levels and 
(^—\^Ko+i Q levels. Inversion accompanied by rotation over tt around one of the 
principle axes corresponds to transposition of the two identical nuclei (H for NH2D, or 
D for ND2H). The allowed levels are determined by the total spin of the identical nuclei. 
Because of that a purely inversion transition is strongly suppressed as an ortho •H- para 
one. 

These selection rules lead to two types of allowed transitions. For transitions with 
= the allowed combination is s •<-)■ a and for transitions with = ±1 the 

allowed combinations are s s and a •<-> a. In other words, transitions with /S.Kc 7^ 
are pure rotational, while transitions with = are mixed. Consequently, all 

transitions with /S.Kc 7^ have — 1, while transitions with /S.Kc — according to 
(10) have sensitivities: 

QM = 1 ± (Qi,,,^ - 1) . (11) 

UJ 

For NH2D molecule the inversion frequency is 12.2 GHz and Qi■av,^I. = 4.7, so 

^ , , 45GHz , , 

Q^{u) = 1 ± . (12) 

UJ 

For ND2H molecule the inversion frequency is a;inv = 5.11 GHz and Qinv.^x = 5.1. Thus, 

^ , , 21GHz , , 

Q^{uj) = l± . (13) 

OJ 

Results for several strong transitions between lowest rotational levels of these molecules 
are given in tables 2 and 3. 
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Table 2. Sensitivity coefBcients for several lower AKc = microwave transitions 
in NH2D. The lower state energy E is given in the third column. 





transition 




E (cm^M 


7/ (GHz) 




li,o(a) 


Oo,o(s) 


para 


0.00 


494.5 


1.09 


li.oGs) 


Oo,o(a) 


ortho 


0.41 


470.3 


0.90 




lo,i(s) 


para 


11.10 


110.2 


1.41 


li,i(s) lo,i(«) 


ortho 


11.51 


85.9 


0.48 


2i,2(a) 


2o,2(s) 


para 


32.78 


74.2 


1.61 


2i,2(s) 


-> 2o,2(a) 


ortho 


33.19 


50.0 


0.10 


22,i(a) 


^2i,i(s) 


ortho 


40.01 


305.7 


1.15 


22,i(s) ^2i,i(a) 


para 


40.41 


282.1 


0.84 



Table 3. Sensitivity coefficients for /S.Kc = transitions in ND2H. 





transition 




E (cm-i) 


V (GHz) 




li,o(a) 


Oo.o(s) 


ortho 


0.00 


388.7 


1.05 


li,o(s) 


-> Oo,o(a) 


para 


0.17 


378.5 


0.94 




lo,i(s) 


para 


9.10 


67.8 


1.31 




U).i{a) 


ortho 


9.27 


57.7 


0.64 


2i,2(a) 


2o,2(s) 


ortho 


26.66 


38.7 


1.54 


2i,2(s) 2o,2(a) 


para 


26.83 


28.6 


0.27 


22,1 (a) 


^2i,i(s) 


ortho 


32.59 


193.0 


1.11 


22,i(s) ^2i,i(a) 


para 


32.76 


183.2 


0.89 



We see that differences between coefficients for ND2H are somewfiat smaller, 
than for NH2D. In combination with significantly smaller abundance, this makes ND2H 
less attractive candidate for the search of //-variation. Table 2 shows that coefficients 
Q/i for AKc = transitions in NH2D molecule with J < 2 vary from 0.1 to 1.6. Thus, 
the largest for deuterated ammonia is almost three times smaller than for the 
inversion line in NH3. However, in NH3 there arc no close lines with different sensitivities 
and we are forced to use rotational lines of other molecules as a reference. For deuterated 
ammonia there are four pairs of lines for ortho and para molecules which have close 
frequencies and excitation temperatures, but significantly different sensitivities. For two 
of these pairs AQ^^ ~ 1. In addition, there are several transitions with AKc — il) which 
all have — 1. This gives very characteristic pattern of frequency variation caused 
by variation of /j, and can be used to suppress various systematic effects. Consequently, 
reliability of the results can be significantly improved. We conclude that microwave 
spectra of deuterated ammonia can be used as important supplement to the ammonia 
method used previously in [3, 4, 5, 6, 7, 8]. 



Sensitivity of deuterated ammonia to variation of electron-to-proton mass ratio 



7 



Acknowledgments 

This research is partly supported by RFBR grants 08-02-00460, 09-02-00352, and 09- 
02-12223. 

References 

[1] Flambaum V V and Kozlov M G 2009 Cold Molecules: Theory, Experiment, Applications (CRC 
Press) chap Variation of Fundamental Constants as Revealed by Molecules: Astrophysical 
Observations and Laboratory Experiments, p 597 {Preprint arXiv: 0711.4536) 

[2] van Veldhoven J, Kiipper J, Bethlem H L, Sartakov B, van Roij A J A and Meijer G 2004 Eur. 
Phys. J. D 31 337 

[3] Flambaum V V and Kozlov M G 2007 Phys. Rev. Lett. 98 240801 {Preprint arXiv : 0704 .2301) 
[4] Murphy M T, Flambaum V V, MuUer S and Henkel C 2008 Science 320 1611 {Preprint 
arXiv: 0806. 3081) 

[5] Henkel C, Menten K M, Murphy M T, Jethava N, Flambaum V V, Braatz J A, Muller S, Ott J 

and Mao R Q 2009 Astron. Astrophys. 500 725 {Preprint arXiv: 0904.3081) 
[6] Levshakov S A, Molaro P and Kozlov M G 2008 On spatial variations of the electron-to-proton 

mass ratio in the Milky Way arXiv: 0808 . 0583 
[7] Molaro P, Levshakov S A and Kozlov M G 2009 Stringent bounds to spatial variations of the 

electron-to-proton mass ratio in the milky way arXiv: 0907. 1192 
[8] Levshakov S A, Molaro P, Lapinov A V, Reimcrs D, Hcnkcl C and Sakai T 2009 Search for 

chameleon-like scalar fields with the ammonia method. Submitted to Astron. Astrophys. arXiv : 

0911.3732 

[9] Burrage C, Davis A C and Shaw D J 2009 Phys. Rev. D 79 044028 {Preprint arXiv: 0809. 1763) 
[10] Kanekar N, CariUi C L, Langston G I et al. 2005 Phys. Rev. Lett. 95 261301 
[11] Levshakov S A, Reimers D, Kozlov M G, Porsev S G and Molaro P 2008 Astron. Astrophys. 479 

719 (Prepri77i arXiv: 0712.2890) 
[12] Kanekar N and Chcngalur J N 2004 Mon. Not. R. Astron. Soc. 350 L17 {Preprint 

arXiv : astro-ph/0310765) 
[13] Kozlov M G 2009 Phys. Rev. A 80 022118 {Preprint arXiv: 0905. 1714) 

[14] Olbcrg M, Bcster M, Rau G, Pauls T, Winnewisser G, Johansson L E and Hjalmarson A 1985 

Astron. Astrophys. 142 LI 
[15] Saito S, Ozeki H, Ohishi M and Yamamoto S 2000 Astrophys. J. 535 227 

[16] Tine S, Roueff E, Falgarone E, Gerin M and Pineau des Forets G 2000 Astron. Astrophys. 356 
1039 

[17] Roueff E, Tine S, Coudert L H, Pineau des Forets G, Falgarone E and Gerin M 2000 

Astron. Astrophys. 354 L63 
[18] Loinard L, Castets A, Ceccarelh C, Caux E and Tielens A G G M 2001 Astrophys. J. 552 L163 
[19] Roueff E, Lis D C, van der Tak F F S, Gerin M and Goldsmith P F 2005 Astron. Astrophys. 438 

585 

[20] Gerin M, Lis D C, Philipp S, Giisten R, Roueff E and Reveret V 2006 Astron. Astrophys. 454 L63 
[21] Lis D C, Gerin M, Roueff E, Vastel C and Phillips T G 2006 Astrophys. J. 636 916 
[22] Lis D C, Gerin M, Roueff E, Phillips T G and Poelman D R 2008 Astrophys. Space Sci. 313 77 
[23] Griest K, Whitmore J B, Wolfe A M, Prochaska J X, Howk J C and Marcy G W 2009 Wavelength 

accuracy of the Keck HIRES spectrograph and measuring changes in the fine structure constant 

arXiv: 0904.4725 

[24] Landau L D and Lifshitz E M 1977 Quantum mechanics 3rd ed (Oxford: Pergamon) 
[25] Swalen J D and Ibers J A 1962 J. Chem. Phys. 36 1914 



